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Abstract 

This paper describes a methodology for performing a temperature retrieval in the Martian 
atmosphere in the 50-90 km altitude range using spectrally integrated 15 pm C0 2 limb emissions 
measured by the Thermal Emission Spectrometer (TES), the thermal infrared spectrometer on 
board the Mars Global Surveyor (MGS). We demonstrate that temperature retrievals from limb 
observations in the 75-90 km altitude range require accounting for the non-local thermodynamic 
equilibrium (non-LTE) populations of the C0 2 (v 2 ) vibrational levels. Using the methodology 
described in the paper, we have retrieved -1200 individual temperature profiles from MGS TES 
limb observations in the altitude range between 60 and 90 km. Our dataset of retrieved temperature 
profiles is available for download in supplemental materials of this paper. The temperature retrieval 
uncertainties are mainly caused by radiance noise, and are estimated to be about 2 K at 60 km and 
below, 4 K at 70 km, 7 K at 80 km, 10 K at 85 km, and 20 K at 90 km. We compare the retrieved 
profiles to Mars Climate Database temperature profiles and find good qualitative agreement. 
Quantitatively, our retrieved profiles are in general warmer and demonstrate strong variability with 
the following values for bias and standard deviations (in brackets) compared to the Martian Year 24 
dataset of the Mars Climate Database: 6 (±20) K at 60 km, 7.5 (±25) K at 65 km, 9 (±27) K at 
70 km, 9.5 (±27) K at 75 km, 10 (±28) K at 80 km, 11 (±29) K at 85 km, and 11.5 (±31) K at 
90 km. Possible reasons for the positive temperature bias are discussed. 
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1. Introduction 


The TES instrument on the MGS orbiter mapped the atmosphere of Mars for more than three 
Martian years allowing an unprecedented survey of the present Martian climate. The TES 
observations have been used to retrieve temperature profiles (Conrath et al. 2000), aerosol optical 
depth and physical properties (e.g. Pearl et al. 2001; Smith et al. 2001; Smith 2004; Clancy et al. 
2003; Wolff and Clancy 2003), and water vapor abundance (Smith 2002; 2004). Most of the 
previously reported temperature retrievals from TES cover the atmosphere up to the -0.01 mbar 
level, which corresponds to -60 km altitude. However, in the limb observation mode measurements 
were performed at higher altitudes, up to -120 km altitude. The upper limit for the currently 
available TES retrievals was set in accordance with the signal-to-noise ratio of the observations 
used in the constrained linear inversion algorithm applied for temperature retrievals (Conrath et al., 
2000). In this work we used a stable forward-fitting temperature retrieval algorithm and applied it 
to a spectrally integrated C0 2 15 pm radiance to extend the upper altitude of temperature retrievals 
to 90 km and to provide the community with high-altitude temperature data for Martian Years 
(MY) 24 and 25. Temperature soundings in the middle atmosphere of Mars are important since 
they can provide important insight into the dynamics of the region that is linked with convective 
instabilities and gravity wave propagation and dissipation (e.g. Heavens et al., 2010 and references 
therein) and can provide additional constraints on tidal and gravity wave drag needed for modeling 
work (e.g. Forget et al., 1999; Forbes and Miyahara, 2006; Hartogh et al., 2007). 

The structure of this paper is as follows: Section 2 describes the TES instrument. In Section 3 we 
discuss the details of the observed C0 2 15 pm band radiance including the breakdown of local 
thermodynamic equilibrium (LTE) at high altitudes and the effects of the upper atmosphere of Mars 
on simulated limb radiances. In Section 3 we also introduce the selection procedure we use on the 
observations and we define the integrated radiance that will be used for the retrievals. Section 4 
describes the non-LTE retrieval code ALI-ARMS (see the Appendix for abbreviations not 
explained in the text), the temperature retrieval algorithm, and the retrieval uncertainties. In 
Section 5 we present the retrieved temperature profiles and compare them to averaged temperature 
profiles from the Mars Climate Database dataset. A discussion of similarities and differences is also 
included into this section. Finally, we present our conclusions in Section 6. 

2. TES instrument 

MGS began mapping operations from Mars orbit on 1 March 1999 (MY 24, Ls = 104°) and 
provided nearly continuous monitoring of conditions in the Martian atmosphere for over three full 
Martian years. The MGS mission ended on 2 November 2006 when the contact with the spacecraft 


was lost. The TES instrument (Christensen et al. 2001) systematically measured and monitored the 
Martian surface and atmosphere throughout all phases of the mission. The imminent failure of a 
neon bulb required for calibration of the spectrometer forced spectrometer observations to end on 
31 August 2004 (MY 27; Ls=81°), although observations continued using the TES bolometers. The 
near-polar solar-locked mapping orbit allowed 12 orbits per day of data at roughly 2:00 AM and 
2:00 PM local time. A majority of observations were taken in the nadir geometry, but limb 
observations of the Martian atmosphere were included roughly every 10 degrees of latitude 
throughout the orbit (both day and night) throughout the mission. Each 8.3 mrad detector field-of- 
view had a projected size of about 13 km at the limb, allowing vertical resolution of just over a 
scale height. In a limb-geometry sequence, the 3x2 array of TES detectors scanned the atmosphere 
in overlapping steps from below the surface to about 120 km tangent height above the surface 
providing vertical sampling of about 1 to 5 km (see for more details in (Christensen et al., 2001)). 
The total number of limb geometry spectra obtained during the mission is about 750,000. In this 
work we show that a limited number of limb geometry spectra taken during the first Martian year 
have signal to noise ratio (SNR) that allows an extension of temperature retrievals up to 90 km 
altitude. The selection criteria for these spectra are described in Section 3.2. 

3. C0 2 15 pm radiance measured by MGS-TES 

In the gas phase, carbon dioxide (C0 2 ) molecular vibrations involve combinations of symmetric 
stretch (Vi), covalent bond bending (v 2 ), and asymmetric stretch (v 3 ) modes. The diagram in Fig. 1 
shows the ground state and various excited vibrational levels of the C0 2 molecule up to an energy 
of 7445 cm -1 . The levels are marked in accordance with the Herzberg notation V|V 2 V. where Vi, v 2 , 
and v 3 denote the number of the corresponding vibrational quanta and the / symbol refers to the 
orbital quantum number. The 15 pm radiance measured by the TES instrument and used for 
temperature retrievals in this study arises from the optical transitions from vibrationally excited 
states with Av 2 =1, where Av 2 denotes the change in the v 2 vibrational quanta number. 


3.1. LTE and non-LTE conditions in C0 2 vibrational levels in the Martian atmosphere 

The interpretation of TES 15 pm limb radiance profiles requires information on the populations of 
the corresponding C0 2 vibrational levels at the altitudes of the limb observations. In the lower 
atmosphere the frequency of inelastic molecular collisions is sufficiently high that these collisions 
overwhelm other population/depopulation mechanisms of the molecular vibrational levels. This 
leads to a local thermodynamic equilibrium, and the populations follow the Boltzmann distribution 
governed by the local kinetic temperature T ki „. In the upper atmosphere of Mars, where the 
frequency of inelastic collisions is much lower than that at lower altitudes, other processes also 
influence the population of C0 2 vibrational levels. These include: a) the direct absorption of solar 
radiance by the C0 2 vibrational-rotational bands in the 2.0, 2.7, and 4.3 pm spectral regions; 

b) absorption of the 15 pm radiance coming from the warmer and denser lower atmosphere; 

c) vibrational-translational (V-T) energy exchanges by collisions with molecules and atoms of 
other atmospheric constituents; and d) collisional vibrational-vibrational (V-V) energy exchange 
with other molecules. As a result, LTE no longer applies in this altitude region and the populations 
must be found by solving the self-consistent system of kinetic and radiative transfer equations, 
which express the balance relations between various excitation/de-excitation processes described 
above. The most important rate coefficient that defines the population of C0 2 (v 2 ) vibrational levels 
in the upper atmosphere of Mars is kvr{C0 2 -0}, the quenching rate coefficient for V-T energy 
exchange process at C0 2 -0 collisions (Ch. 10.3 in (Lopez-Puertas and Taylor, 2001)). Despite the 
efforts of several different research groups, there is still a discrepancy between the laboratory 
measurements of this rate coefficient that provide low values of about 1.3-1. 5 x 10~ 12 cm 3 s _1 
(Pollock et al., 1993; Khvorostovskaya et al., 2002; Castle et al., 2006) and atmospheric retrievals 
that require kvr{C0 2 -0} = 6.0 x 10 12 cm’s 1 for an adequate interpretation of 15 pm C0 2 radiance 
(Sharma and Wintersteiner, 1990; Ratkowski et al., 1994; Gusev et al., 2006). The discussion of 
this problem is beyond the scope of this work, but its effect is minimized in this work because the 
C0 2 (v 2 ) vibrational level populations are in LTE up to ~80-90 km, and we utilize a “buffer zone” 
(see the end of Sect. 3.1 and Sect. 4) in the 90-100 km altitude range that compensates for the 
uncertainties in the upper part of the atmosphere including uncertainty in the k V T{C0 2 -0} rate 
coefficient. The value kvr{C0 2 -0} = 6.0 x 10 l2 cm’s 1 was used in all calculations performed for 
this work. 

In this retrieval we used the non-LTE model developed in our group (Gusev, 2003; Kutepov et al., 
2006). This model includes five C0 2 isotopes, 60 vibrational levels, and 20,000 optical transitions 


from the HITRAN 2004 database (Rothman et al. 2005). It can easily be extended to include more 
levels and more transitions; however, this was not needed for the purposes of our work (see the 
discussion of the levels needed to be included into the task in Ch. 10.3 of (Lopez-Puertas and 
Taylor, 2001)). The calculations and retrievals were performed using the ALI-ARMS research non- 
LTE code (see Section 4 below). 

The calculations show that the main contributors to 15 pm radiance are the fundamental and first 
hot 15 pm C0 2 bands and that the corresponding 01 J 0 and 02 2 0 vibrational levels are in LTE below 
~85 km altitude. Even though the area of our primary interest (60-90 km) doesn’t demonstrate 
strong non-LTE features, performing non-LTE calculations is still required to take into account the 
contributions of the upper layers in the limb observation geometry. This is demonstrated in Fig. 2, 
where a reference profile of C0 2 15 pm radiance is compared with the C0 2 radiance profiles 
calculated using different approaches. The reference profile of non-LTE 15 pm radiance was 
estimated for the tropical atmosphere (lat = 0°), L s =180° model extracted from the Mars Climate 
Database (MCD). We used pressure/temperature, C0 2 , and O VMR profiles covering the 0-200 km 
altitude range. In this and other exercises we used the recommended MY 24 dust and average solar 
flux scenario dataset available for download at 

http://johnson.lmd.jussieu.fr:8080/las/servlets/dataset. Here and below we will refer to this dataset 
as simply the “MCD”). To obtain this profile we first calculated the non-LTE populations of the 
C0 2 vibrational levels at all heights and then calculated the limb radiance for each spectral line in 
the 560-770 cm 1 spectral interval and at each tangent height. After that, the radiance was spectrally 
integrated and the resulting integrated 15 pm radiance formed the reference profile. To demonstrate 
the importance of the accurate calculation of C0 2 (v 2 ) populations in the upper atmosphere for the 
limb radiance simulations below, we simulated limb radiance using LTE populations of vibrational 
levels (curve with circles in Fig. 2). As one can see, this significantly overestimates the radiance. In 
the course of the retrieval this approach would lead to underestimating the temperature down to 
50 km altitude. One of the ways of dealing with this problem is by cutting out the upper atmosphere 
in the calculations to avoid the unrealistic radiance coming from the upper layers. Figure 2 shows 
the results of LTE and non-LTE runs when the atmosphere was cut out just above 100 km. As one 
can see, in these cases the radiance is underestimated above ~80 km altitude, which would lead to 
overestimating temperatures in the retrieval. Moving the upper limit of the atmosphere in the 
calculations, one can compensate for the non-LTE correction. However, this will introduce an 
uncontrolled uncertainty since the upper limit will depend on the atmospheric profile and individual 
non-LTE calculations will be required to find the corresponding altitude. On the other hand, using 
the full non-LTE model in combination with averaged atmospheric profiles in the upper 


atmosphere helps to eliminate this problem. In addition (see Section 4), in our temperature 
retrievals we used the 90-100 km altitude interval as a “buffer zone” to compensate for possible 
deviations of the upper atmosphere profiles from the real ones. Temperature retrievals in the buffer 
zone are not included to the output product and are presented in data files only for control purposes: 
the stronger the gradient in the buffer zone, the larger the discrepancy of the simulated and 
measured profiles above 100 km altitude due to the uncertainties in this area. 

3.2. MGS-TES spectra selection and processing 

In this work we used spectrally integrated 15 pm C0 2 band radiance in the retrievals to increase 
SNR and to extend the upper limit of the retrieved temperature profile. In a sense, using the 
integrated radiance converts the spectrometer to a broadband radiometer like, for example, the 
SABER instrument (Russell et al., 1999) that was successfully used for temperature retrievals in the 
Earth’s atmosphere up to ~110 km altitude. Spectra selection for integration and processing was 
based on two criteria: a limb sequence should a) provide a sufficient number of measured points 
above 60 km altitude and b) SNR for the integrated signal should be greater than unity up to at least 
85 km. The calculations performed in this work show that both during the day and night in the 
Martian atmosphere, the 15 pm C0 2 band radiance is confined within the 560-770 cm' 1 spectral 
range and all spectral line intensities in the 510-560 cm" 1 and 770-830 cm' 1 spectral windows are 
less than 1% of the line intensities in the maximum of the 15 pm manifold. Correspondingly, we 
have selected the former interval as the “information content” window and the two latter as “noise 
and bias” windows. To remove biases we fit for a second order polynomial using the spectral points 
in the “noise and bias” windows, and then we subtract the corresponding values from the 
“information content” spectral intensities. Examples of the resulting spectra in the 61-89 km 
altitude interval are shown in Fig. 3a-f. As one can see, the suggested procedure efficiently 
removes any bias at all considered altitudes. The integrated values of the 15-pm radiance in the 
560-770 cm' 1 spectral range are shown in Fig. 3g. This panel also contains the NER value 
estimated from observed signal values in the “noise and bias” windows. As one can see in Fig. 3g, 
the NER becomes comparable to the integrated signal at ~86 km altitude. The number of limb 
sequences for which the NER allows retrievals up to 90 km altitude is small compared with the 
total number of limb measurements. There are two reasons for that. First, the altitude coverage of 
limb measurements differs and not all limb sequences provide valid data points up to 90 km 
altitude. Second, the increasing noise in TES observations as the mission went on limits the limb 
scans with “good SNR” up to 90 km to the first two Martian years of the mission. Selecting high 
SNR scans may introduce a positive bias to the average retrieved temperature profile as we will 
show below in Sect. 5.2. However, this does not produce any errors in individual temperature 


profiles that we have used for analysis. Overall, we have identified f410 limb sequences that satisfy 
the SNR and altitude coverage criteria (see Fig. 4). All of these observations were processed using 
the retrieval procedure described in the next section. 

4. ALI-ARMS research code and retrieval algorithm 

All of the calculations performed in this work were made using the ALI-ARMS computer code (see 
Kutepov et al., 1998; Gusev, 2003; Gusev and Kutepov, 2003; and references therein) that solves 
the multi-level problem using the Accelerated Lambda Iteration (ALI) technique developed for 
calculating non-LTE populations of atomic and ionic levels in stellar atmospheres (Rybicki and 
Hummer, 1991). The code iteratively solves a set of statistical equilibrium equations and the 
radiative transfer equations. The algorithm efficiency is ensured by the ALI technique, which 
avoids the expensive radiative transfer calculations for the photons trapped in the optically thick 
cores of spectral lines. The ALI-ARMS model was successfully applied by Kaufmann et al. (2002; 
2003) and Gusev et al. (2006) to the non-LTE diagnostics of spectral observations of the Earth’s 
limb from the CRISTA instrument (Offermann et al., 1999; Grossmann et al., 2002). Kutepov et al. 
(2006) used this model to validate temperature retrievals from the 15 pm C0 2 emissions measured 
by the SABER instrument onboard the terrestrial TIMED satellite. Feofilov et al. (2009) used this 
code for the H 2 0 non-LTE model validation and H 2 0 VMR retrievals from the 6.3 pm H 2 0 
emissions measured by SABER. 

The retrieval method implemented in the ALI-ARMS code is similar to the iterative onion-peel 
technique using the relaxation method described in Gordley and Russell (1981). The process starts 
with an initial guess of a temperature profile combined with a fixed atmospheric model (pressure, 
VMRs of C0 2 and atomic oxygen retrieved from the MCD). The non-LTE populations are 
calculated and used for monochromatic limb radiance calculations for each limb path that is 
integrated and convolved with the MGS-TES field of view. The resulting simulated radiance I is 
compared to the measured radiance at each tangent height, and the temperature profile is iterated 
using the following relaxation scheme: £‘ +1 = % 1 + ( I meas -!')/( dl /d%), where C, 1+1 and £ ' are the 
temperature values at the z+i-th and z-th iterations, respectively, I meas is the integrated limb radiance 
measured by TES (see Section 3.2 above), 1 1 is the simulated limb radiance at the z-th iteration, and 
( dl /dc, ) is the numerically calculated derivative of the radiance produced by the forward model 
with respect to After all limb paths have converged, a new temperature profile is produced, the 
hydrostatic equilibrium equation is solved for, the new non-LTE populations of C0 2 molecular 
levels are calculated, and the radiance is simulated again. The iterations are repeated until the 


differences between the simulated and measured radiances become equal to or smaller than the 
radiance noise in the observed integrated limb radiance. 

Figure 5 shows a self-consistency check of the retrieval procedure performed for the model 
atmosphere described above. First, the non-LTE task was solved for the reference temperature 
profile T re ^z), and the limb radiance was calculated using the non-LTE populations of the C0 2 
levels, observation geometry, and TES field of view. Second, the resulting radiance profile was 
used as the “measured” profile / mras (z). Two tests were then conducted. The initial guess 
temperature profile T(z) was set equal to T re f(z ) - 30 K (case I) and T re f(z) + 30 K (case II) in the 
60 - 90 km altitude range. Below 50 km and above 100 km the original values of the T re f(z) profile 
were used. In the 50-60 km and 90-100 km altitude ranges the temperature profiles were smoothly 
sewn together. The retrieval procedure was run for both cases. Figure 5a shows that the retrieved 
temperature profile rapidly converges to the reference profile in the course of a few iterations and 
that the result only slightly depends on the initial profile. This dependence is due to field of view 
averaging effects. To minimize the temperature uncertainty caused by these effects the retrieval 
procedure should comply with the following requirements: a) the initial guess T(z) profile should 
not be unrealistically high or low. We recommend starting the retrieval from the climatological 
mean T(z). b) The iterations should be stopped as soon as the convergence criterion is reached. 
Further fitting of radiance may lead to the appearance of false structures in the retrieved 
temperature profiles (Houghton et al, 1984, p.148). Figure 5b demonstrates the differences between 
the simulated and reference radiances after each iteration. One can see that the converged radiance 
profile reproduces the reference profile well at all points in the 60 - 90 km altitude range. We also 
performed runs similar to those described above but with noise added to I meas ( z) and for different 
seasons and geographic locations. Noise levels were estimated from real data. This study provided 
the following temperature retrieval uncertainties: 2 K at 60 km and below, 4 K at 70 km, 7 K at 
80 km, and 10 K at 85 km altitude. 

5. Temperature retrievals from TES measurements 

5.1. Input data and retrieval procedure 

All the retrievals performed in this study utilized the retrieval algorithm described in Section 4. For 
each of the 1410 integrated TES radiance profiles (see Section 3.2) the initial guess atmospheric 
profile was assembled from pressure/temperature profiles in 10-60 km altitude range retrieved by 
the TES team using nadir geometry observations (Conrath et al. 2000; Smith 2004), and from 
pressure/temperature, C0 2 and O VMRs taken from the MCD in the 60-200 km interval. For each 
integrated radiance profile the retrieval procedure was run for 20 retrieval iterations. The retrieval 
was performed over the 50-100 km range with the 90-100 km interval used as a “buffer zone” to 


smooth out the effects of limb radiance discrepancies in the upper atmosphere (see Section 3.1). 
The 50-60 km interval was included in the retrieval to allow for temperature profile adjustment in 
this area due to changes in the upper part. For each of the 20 iterations the root mean square (r.m.s.) 
of the radiance deviation ( I meas - Laic ) was calculated in the 60-90 km altitude range. For the 
majority of the profiles the minimum r.m.s. radiance deviation was reached between the 7th and 
10th iteration. If convergence was not reached or the deviation of radiances at any altitude of the 
retrieval interval range was greater than 10% then the profile was not included to the output data. 
Approximately 14% of the profiles did not pass this test. We present the remaining 1214 
temperature profiles in this work. 

5.2. Retrieved temperature profiles 

The complete set of individual temperature profile retrievals is available for download as an archive 
at http://cua-nasa-gsfc.info/feofilov/TES_HighAlt.zip. For a description of the dataset see the file 
“readme.txt” inside the archive. In this work we present some examples of individual retrievals and 
compare the averaged profiles with those from the MCD (using the standard MY 24 dust and 
average flux scenario). 

Three examples of individual temperature profiles retrieved from integrated TES radiance are 
shown in Fig. 6a-c. As expected, below ~55 km they match the initial temperature profile. Above 
60 km, the retrieved temperature profiles differ from the model and, in general, individual 
temperature retrievals show much higher variability than the model profiles. This may be explained 
by gravity waves in the Martian atmosphere, which are not reproduced in the model. Lack of 
temperature measurements in the Martian atmosphere for the period considered in the study 
prevents a profile-by-profile comparison with other instruments. Instead, we compare the zonally 
and seasonally (with a 30° step over L s ) averaged temperature profiles from TES retrievals with 
corresponding profiles from the MCD. We present 12 averaged temperature profiles in Figure 7. 
The profiles are marked in the figure caption with L s counted from the year the mission has started 
(MY 24) so that L s =390° corresponds to MY 25, L s =30°. Each of the rows of panels in Fig. 7 
contains three temperature distributions: high latitudes in the Southern Hemisphere, equatorial or 
mid-latitude, and high latitudes in Northern Hemisphere. Each displayed profile represents an 
average over ~10-15 individual retrieved profiles. As one can see, there is qualitative agreement 
between the retrieved temperature profiles and the MCD. However, absolute values differ and in 
some cases the difference may exceed 30 K. In general, our TES retrievals are warmer in the 60- 
90 km altitude range than the MCD temperature profiles. An analysis of our entire retrieval set in 
comparison to the MCD profiles gives the following values for bias and standard deviations (in 
brackets): 6 (±20) K at 60 km, 7.5 (±25) K at 65 km, 9 (±27) K at 70 km, 9.5 (±27) K at 75 km, 


10 (±28) K at 80 km, 11 (±29) K at 85 km, and 11.5 (±31) K at 90 km. The temperature variability 
responsible for the standard deviations is expected since our TES retrievals represent individual 
profiles. The increase of standard deviation with altitude can be explained by the growth of gravity 
wave amplitude (Fritts et al., 2006). The positive temperature bias can be explained by a number of 
reasons. First, there might be deviations of the real atmospheric conditions for a given period from 
that modeled by the GCM. It is known that there were at least two major dust storm events during 
the considered time period (Smith, 2004): one at L s =225° and one at L s =260° that might be the 
cause for large temperature deviation shown in Fig. 7b. Second, the positive temperature bias may 
be caused by offsets in the TES retrievals. One cannot assign these offsets to uncertainties of the 
non-LTE model since those only affect the upper part of the retrieval domain. One should also rule 
out the retrieval algorithm issues since, as the self-consistency study shows (see Section 4), it can 
introduce only minor disturbances to the temperature profile. The offset may be caused by radiance 
calibration issues like the vertical field of view wing shape or by the spectra selection procedure. A 
radiance leak from below would be barely noticeable if the vertical radiance gradients are low. 
However, if an area with low temperatures and pressures is observed then the contributions of 
lower layers become important (Rezac, 2011). Positive bias in our TES retrievals can also be linked 
with the limb sequence selection criteria described in Section 3.2. Using a selecting procedure that 
discards the scans with low SNR may also filter out the low temperature cases that are 
characterized by low radiance values. 

6. Conclusion 

We have described a non-LTE model and algorithm that is applied to the temperature retrieval 
using the 15 pm limb geometry emissions measured by the TES instrument. We used integrated 
limb radiance in the 560-770 cm' 1 spectral range to increase the SNR and to extend the retrievals 
up to ~90 km altitude. We applied a forward fitting retrieval algorithm and showed that it provides 
a stable solution that does not depend on the initial guess profile. The combined uncertainty of the 
temperature retrieval due to noise and algorithm implementation is estimated to be 2 K at 60 km 
and below, 4 K at 70 km, 7 K at 80 km, 10 K at 85 km, and 20 K at 90 km altitude. Performing the 
retrievals from the integrated TES radiance produces temperature profiles similar to those measured 
by other instruments and predicted by models. Overall, ~1200 individual limb sequences have been 
selected, processed, and provided for download by the community. The temperature profiles 
retrieved from TES at 60 - 90 km altitude were compared to corresponding profiles from the Mars 
Climate Database. In general, good qualitative agreement between the two datasets was observed. 
A positive bias of TES retrieved temperatures with respect to MCD temperatures was found that 
varies from 6 K at 60 km to 11.5 K at 90 km. We explain this bias both by natural atmospheric 


variability and by the filtering out of some low temperature profiles in our signal selection routine. 
The standard deviation of our retrieved temperature profiles compared to the MCD varies from 
±20 K at 60 km to ±31 K at 90 km, which is explained by gravity wave growth with altitude as well 
as by larger retrieval uncertainties at higher altitudes. In summary, this study shows how to extend 
the limits of the already highly capable TES instrument, and validates the general concepts built in 
to Martian GCMs. 
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Appendix: Abbreviations used in text 

ALI-ARMS: Accelerated Lambda Iterations for Atmospheric Radiation and Molecular Spectra 
CRISTA: CRyogenic Infrared Spectrometers and Telescopes for the Atmosphere 
GCM: General Circulation Model 

HITRAN: High-resolution TRANsmission molecular absorption database 
LTE: Local Thermodynamic Equilibrium 
MGS: Mars Global Surveyor 

MY24: Martian Year 24, Mars Climate Database scenario characterized by low dust conditions 
NER: Noise Equivalent Radiance 

SABER: Sounding of the Atmosphere using Broadband Emission Radiometry 

SNR: Signal to Noise Ratio 

TES: Thermal Emission Spectrometer 

TIMED: Thermosphere Ionosphere Mesosphere Energetics and Dynamics 
VMR: Volume Mixing Ratio 
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Fig. 1. Scheme of vibrational level and major energy exchange processes for C0 2 and its isotopes. 
Solid lines: radiative transitions, dashed lines: V-V energy transfer from solar pumped levels. V-T 
processes are not shown. FB: fundamental band; FH: first hot band; SFi: second hot band. 

Fig. 2. Changes in the integrated f 5 pm radiance using the L s = 180°, lat = 0° Mars Climate 
Database atmospheric profile. Test runs were performed for LTE conditions in the atmosphere set 
up to 200 km (circles), non-LTE conditions in the atmosphere cut off at 100 km (triangles), and 
LTE conditions in the atmosphere cut off at 100 km (squares). 

Fig. 5. Self-consistent retrievals of the temperature profile T re f(z) in the 60 - 90 km altitude range: 
a) iterations starting with two different initial guess temperature profiles: T(z) = T, af (z) - 30 K (case 
I) and T(z) = T ref (z) + 30 K (case II); b) percent difference between the reference and calculated 
radiance profiles in the course of iterations for cases I and II. 

Fig. 3. MGS TES spectra over 530-770 cm 1 at: a) 88.6 km; b) 82.6 km; c) 75.6 km; d) 71.4 km; 
e) 68.2 km; f) 60.8 km. Panel g) shows the 15 pm radiance integrated over the 560-772 cm ' 
spectral range. Noise equivalent radiance (NER) was estimated from the points outside the 
integration interval: 508-560 cm ' and 772-825 cm 4 . 

Fig. 4. Seasonal and latitudinal coverage of the MGS TES limb sequences satisfying the selection 
criteria. 

Fig. 6. Temperature retrieval examples: a) L s =110°, lat = -62°; b) L s =135°, lat = -6°; c) L s =180°, 
lat = 59°. Curves with circles: initial temperature profile (MGS-TES below 50 km; MCD above 
50 km). Curves with triangles: temperature retrievals performed in this work. Dashed lines: error 
bars for temperature retrievals. 

Fig. 7. Zonally (±2.5°) and seasonally (±15°) averaged temperature retrievals from MGS-TES 
integrated 15 pm C0 2 radiances. Curves with circles: reference temperature profile (MGS-TES 
below 50 km; MCD above 50 km). Curves with triangles: temperature retrievals performed in this 
work: a) L s =240°, lat = -73°; b) L s =240°, lat = -7° (note the larger temperature difference perhaps 
due to a dust storm); c) L s =240°, lat = 64°; d) L s =270°, lat = -79°; e) L s =270°, lat = 25°; f) L s =270°, 
lat = 64°; g) L s = 360°, lat = -68°; h) L s =360°, lat = 20°; i) L s =360°, lat = 69°; j) L s =390°, lat = - 
88°; k) L s =390°, lat = -24°; 1) L s =390°, lat = 64°. 


